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Resonant proximity effect in the normal metal / insulator / diffusive ferromagnet / insulator / 
s- and d-wave superconductor (N/I/DF/I/S) junctions is studied for various regimes by solving the 
Usadel equation with the generalized boundary conditions. Conductance of the junction and the 
density of states in the DF layer are calculated as a function of the insulating barrier heights at the 

■ interfaces, the magnitudes of the resistance, Thouless energy and the exchange field in DF and the 
' misorientation angle a of a d-wave superconductor. It is shown that the resonant proximity effect 

originating from the exchange field in DF layer strongly modifies the tunneling conductance and 
density of states. We have found that, due to the resonant proximity effect, for s-wave junctions a 
sharp zero bias conductance peak (ZBCP) appears for small Thouless energy, while a broad ZBCP 
appears for large Thouless energy. The magnitude of this ZBCP can exceed the normal state 
conductance in contrast to the case of diffusive normal metal / superconductor junctions. Similar 
structures exist in the density of states in the DF-layer. For d-wave junctions at a = 0, similar 
structures are predicted in the conductance and the density of states. With the increase of the angle 
a, the magnitude of the resonant ZBCP decreases due to the formation of the mid gap Andreev 
resonant states. 



I. INTRODUCTION 



There is a continuously growing interest in the physics of charge and spin transport in ferromagnet / superconductor 
(F/S) junctions. One of the applications of F/S junctions is determination of the spin polarization of the F layer. 
Analyzing signatures of Andreev reflection 1 - in differential conductance by a modified Blonder, Tinkham and Klapwijk 
(BTK) theory^, one can estimate the spin polarization of the F \&ye£££&i2&. This method was generalized and 
applied to ferromagnet / unconventional superconductor junctions 9 . Most of these works are applicable to ballistic 
ferromagnets while understanding of physics in contacts between diffusive ferromagnets (DF) and (both conventional 
and unconventional) superconductors (S) is not complete yet. The model should also properly take into account the 
proximity effect in the DF/S system. 

In DF/S junctions Cooper pairs penetrating into the DF layer from the S layer have nonzero momentum due to the 
exchange field 10 i n i 12 i 13 i 14 i 15 . This property results in many interesting phenomena 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 i 29 i 30 . 
\ One interesting consequence of the oscillations of the pair amplitude is spatially damped oscillating behavior of the 
density of states (DOS) in a ferromagnet predicted theoretically^i* 3 ^ 3 ^* 3 ^ . In a strong ferromagnet the exchange 
field breaks the induced Cooper pairs, while for a weak exchange field the pair amplitude can be enhanced and the 
energy dependent DOS can have a zero-energy peakiSSMiSiSiiS 7 .. Since DOS is a fundamental quantity, this resonant 
proximity effect can influence various transport phenomena. In our recent paper— the DOS peak was studied in two 
regimes of weak and strong proximity effect and the conditions for the appearance of this DOS anomaly were clarified. 
However, its consequence for the junction conductance was not systematically investigated so far. 

It is known that in contacts involving unconventional superconductors the so-called zero-bias conductance peak 
(ZBCP) takes place due to the formation of the midgap Andreev resonant states (MARS)22i22i!ii— An interplay of 
the resonant proximity effect with MARS in DF/d-wave superconductor (DF/D) junctions is an interesting subject 
• i-h , which deserves theoretical study. 

The purpose of the present paper is to formulate theoretical model for the proximity effect in the normal 
metal/DF/s- and d-wave superconductor (N/I/DF/I/S) junctions and to study the influence of the resonant proximity 
effect due to the exchange field on the tunneling conductance and the DOS. A number of physical phenomena may 
coexist in these structures such as impurity scattering, oscillating pair amplitude, phase coherence and MARS. We 
will employ the quasiclassical Usadel equations^ with the Kupriyanov-Lukichev boundary conditions 4 * generalized 
by Nazarov within the circuit theory 4 ^. The generalized boundary conditions are relevant for the actual junctions 
when the barrier transparency is not small. New physical phenomena regarding zero-bias conductance are properly 
described within this approach, e.g., the crossover from a ZBCP to a zero bias conductance dip (ZBCD). The gen- 
eralized boundary conditions were recently applied to the study of contacts of diffusive normal metals (DN) with 
conventional 4 ^ and unconventional superconductors"^. Here we consider the case of N/I/DF/I/S junctions with 
a weak ferromagnet having small exchange field comparable with the superconducting gap. SF contacts with weak 
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ferromagnets were realized in recent experiments with, e.g., CuNi alloys^, Ni doped PdS or magnetic semiconductors. 
Therefore, our results are applicable to these materials and may be observed experimentally. 

The normalized conductance of the N/I/DF/I/S) junction ar{eV) = as(eV)/aN(eV) will be studied as a function 
of the bias voltage V, where 0s(JV) ( e ^0 is the tunneling conductance in the superconducting (normal) state. We will 
consider the influence of various parameters on <7y(eV), such as the height of the interface insulating barriers, the 
resistance Rd, the exchange field h and the Thouless energy Etk in the DF layer. In the case of d-wave superconductor, 
important parameter is the angle between the normal to the interface and the crystal axis of <i-wave superconductor 
a. Throughout the paper we confine ourselves to zero temperature and put ks = Ji = 1- 

The organization of this paper is as follows. In section II, we will provide the detailed derivation of the expression 
for the normalized tunneling conductance. In section III, the results of calculations are presented for various types of 
junctions. In section IV, the summary of the obtained results is given. 



II. FORMULATION 



In this section we introduce the model and the formalism. We consider a junction consisting of normal and 
superconducting reservoirs connected by a quasi-one-dimensional diffusive ferromagnet conductor (DF) with a length 
L much larger than the mean free path. The interface between the DF conductor and the S electrode has a resistance 
Rb while the DF/N interface has a resistance R' b . The positions of the DF/N interface and the DF/S interface are 
denoted as x = and x = L, respectively. We model infinitely narrow insulating barriers by the delta function U{x) — 
H5(x — L)+H'S(x). The resulting transparency of the junctions T m and T' m are given by T m = 4cos 2 </>/(4cos 2 4>+Z 2 ) 
and T' m — 4 cos 2 </)/(4cos 2 </> + Z' ), where Z — 2H/vp and Z' = 2H'/vf are dimensionless constants and </> is the 
injection angle measured from the interface normal to the junction and vp is Fermi velocity. 

We apply the quasiclassical Keldysh formalism in the following calculation of the tunneling conductance. The 
4x4 Green's functions in N, DF and S are denoted by Gq(x), Gi(x) and G 2 {x) respectively where the Keldysh 
component Kq\<x(x) is given by K~i(x) — Ri{x)f i {x) — f i (x)A i {x) with retarded component Ri{x), advanced component 
Ai(x) — —R*(x) using distribution function fi(x)(i = 0,1,2). In the above, Ro(x) is expressed by Ro(x) = T3 
and fo(x) = f w + f 3 f t0 . R 2 {x) is expressed by R 2 {x) = gr 3 + ff 2 with g = e/Ve 2 - A 2 and / = A/VA 2 - e 2 , 
where t 2 and f 3 are the Pauli matrices, and e denotes the quasiparticle energy measured from the Fermi energy 
and f 2 {x) — tanh(e/2T) in thermal equilibrium with temperature T. We put the electrical potential zero in the 
S-electrode. In this case the spatial dependence of G\ (x) in DF is determined by the static Usadel equation^ 3 ., 



D ^ [Gl{x) ^M ]+i[ H,G 1 (x))=0 (1) 



dx dx 

with the diffusion constant D in DF. Here H is given by 

V Ho 

with Hq = (e — (+)/i)t3 for majority(minority) spin where h denotes the exchange field. Note that we assume a weak 
ferromagnet and neglect the difference of Fermi velocity between majority spin and minority spin. The Nazarov's 
generalized boundary condition for G\(x) at the DF/S interface is given in Refs4§4&. The generalized boundary 
condition for G\{x) at the DF/N interface has the form: 

A (6l §) |a=0+ = -^< s y, (2 ) 



B= 2r i UG (0_),G 1 (0 + )] 

4 + T4([G (0_),G 1 (0 + )] + -2)- 

The average over the various angles of injected particles at the interface is defined as 



rn/2 
(,)_ J-tt/2 ' 



/:/ 2 /2 d0T(')(0)cos0 



3 



with B(cj>) — B and T^{(f>) — Tm • The resistance of the interface Rb is given by 



° /!f /2 ^TW(0)cos^' 

Here i?Q^ is Sharvin resistance given by R^- 1 = e 2 fc|S'W/(47r 2 ) in the three-dimensional case. 
The electric current per spin direction is expressed using G\ (x) as 



^(AM^A (3) 

where (Gi(a:) 8G g 1 ^ ) ;i ' denotes the Keldysh component of (Gi(x) 9G q^ )■ In the actual calculation it is convenient 

to use the standard ^-parameterization where function R\{x) is expressed as R\{x) = t^cos9{x) + T2SmO(x). The 
parameter 6{x) is a measure of the proximity effect in DF. 

The distribution function f\{x) is given by f\{x) = fi(x) + f^f t {x) where the component ft{x) determines the 
conductance of the junction we are now concentrating on. From the retarded or advanced component of the Usadel 
equation, the spatial dependence of 6{x) is determined by the following equation 

D^6{x)+2i(e- (+)h)sm[9(x)] =0 (4) 
for majority(minority) spin, while for the Keldysh component we obtain 

£>^[^^cosh 2 ^ m (a ; )]=0. (5) 
At x = 0, since fto is the distribution function in the normal electrode given by 

fto = ~{tanh[(e + eV)/(2T)] - tanh[(e - eV)/(2T)}}. 
Next we focus on the boundary condition at the DF/N interface. Taking the retarded part of Eq. we obtain 



L d0(x) < F >' 

~R d ^x~ U=0+ " 



(6) 



h 



2T' sin0 o 



(2-T; n )+T/ n cos0 o ' 



with 6» = 6»(0+) . 

On the other hand, from the Keldysh part of Eq. (^J) , we obtain 



L ,df t 2 <jq >' (/ t0 -/ t (0 + )) 

-( — )cosh flim (x) U =0+ = ^ , (7) 



with 



T = T£A[ + 2r/„(2 - T^)Real{cosg } 
|(2-T4)+T4cos^o| 2 



A; = (1+ | cos6» | 2 + | sin6l | 2 ). 
Finally, we obtain the following final result for the electric current through the contact 



_j_ f°° ^ e fto 

<ho> L JO cosK 2 9 im (x) ^ </ M >' 
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Then the differential resistance R per one spin projection at zero temperature is given by 

R = ^- + 2 -^[ L *E + ^L_ (9) 



with 

ho = 



T I 2 „A 1 + 2T m (2-T m )A 2 



2 | (2-T m ) + T m [g cos 9 L + fsm6 L ] | 2 ' 
Ax = (1+ | co S 9 L | 2 + | sin0 L | 2 )(| 5 | 2 + | / | 2 +1) 

+4Imag[/<7*]Imag[cos0 i sin0£], (10) 

A 2 = Real{g(cos0 L + cos(9£) + f(sin0 L + sin0£)}. (11) 

This is an extended version of the Volkov-Zaitsev-Klapwijk formula^. For a d-wave junction, the function I^q is 
jiven by the following expression^ 

T T m Cq 
ho 



2 | (2 - T m )(l + g + g_ + /+/_) + T m [cos^( ff+ + <?_) + sm6 L (f+ + /_)] | 2 
Co = T m (l+ | cos0 L | 2 + | sin0 L | 2 )[| 5+ + g_ | 2 + | /+ + /_ | 2 + | 1 + /+/_ + ,g+.g_ | 2 + | /+.g_ - g + f_ | 2 ] 
+2(2 - T m )Real{(l + 5 ;.g* + /;.D[(cos f? L + cos 9* L ) (g+ + gJ) + (sin0 L + sin^)(/+ +/_)]} 
+4T m Imag(cos0 i sin0£)lmag[(/+ + /_)(<£ + <?*)], 



g± — s/y£ 2 — A±, /± = A±/y — e 2 and A± = A cos 2(0 =p a). In the above a, 9 im (x) and 0^ denote the 

angle between the normal to the interface and the crystal axis of d-wave superconductors, the imaginary part of 
0(x) and 0(L-) respectively. Then the total tunneling conductance in the superconducting state as(eV) is given by 
as(eV) = 1 1/R. The local normalized DOS N(s, x) in the DF layer is given by 



N(e,x) = -^Recos6»(x). 



2 

T4 

It is important to note that in the present approach, according to the circuit theory, Rd/R^p can be varied in- 
dependently of Tm\ i.e., independently of Z^'\ Based on this fact, we can choose Rd/R^ 1 and as independent 
parameters. 

In the following section, we will discuss the normalized tunneling conductance <tt(sV) = &s{eV)/ a 'N{eV) where 
tjN(eV) is the tunneling conductance in the normal state given by <7jv(eV) = = l/(-R<2 + Rb + R'b)- 

III. RESULTS 

In this section, we study the influence of the resonant proximity effect on tunneling conductance as well as the DOS 
in the DF region. The resonant proximity effect was discussed in Ref. 3 - and can be characterized as follows. When 
the proximity effect is weak (Rd/Rb <C 1), the resonant condition is given by Rd/Rb ~ ^h/Erh due to the exchange 
splitting of DOS in different spin subbands. When the proximity effect is strong (Rd/Rb 3> 1), the condition is given 
by Eth ~ h and is realized when the length of a ferromagnet is equal to the coherence length £p = y/D/h. We choose 
Rd/Rb = 1 and Rd/Rb = 5 as typical values representing the weak and strong proximity regime, respectively. We fix 
Z' = 3 because this parameter doesn't change the results qualitatively and consider the case of high barrier at the 
N/DF interface, Rd/R' b — 0.1, when the proximity effect is strong. 
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A. Junctions with s-wave superconductors 

We first choose the weak proximity regime and relatively small Thouless energy, EtkI 'A = 0.01. In this case the 
resonant condition is satisfied for h/A — 0.005. In Fig. ^ the tunneling conductance is plotted for Rd/Rb = 1, 
E Th /A = 0.01 and various h/A with (a) Z = 3 and (b) Z = 0. The ZBCP and ZBCD occur due to the proximity 
effect for h — 0. For h/A = 0.005, the resonant ZBCP appears and split into two peaks or dips at eV ~ ±h with 
increasing h/A. The value of the resonant ZBCP exceeds unity. Note that ZBCP due to the conventional proximity 
effect in DN/S junctions is always less than unity^S^fl*Si and therefore is qualitatively different from the resonant 
ZBCP in the DF/S junctions. 

The corresponding normalized DOS of the DF is shown in Fig. 2. Note that in the DN/S junctions, the proximity 
effect is almost independent on Z parameter. We have checked numerically that this also holds for the proximity 
effect in DF/S junctions. Figure displays the DOS for Z = 3, R d /Rb = 1 and E T h/A = 0.01 with (a) h/A = and 
(b) h/A = 0.005 corresponding to the resonant condition. For h = 0, a sharp dip appears at zero energy over the 
whole DF region. For nonzero energy, the DOS is almost unity and spatially independent. For h/A = 0.005 a zero 
energy peak appears in the region of DF near the DF/N interface. This peak is responsible for the large ZBCP shown 
in Fig. ^ Therefore ZBCP in DF/S junctions has different physical origin compared to the one in DN/S junctions. 
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FIG. 1: (color online) Normalized tunneling conductance for s-wave superconductors with Rd/Rb = 1 and Eth/A = 0.01. 

Next we choose the strong proximity regime and relatively small Thouless energy, Eth/A — 0.01. In the present 
case, the resonant ZBCP is expected for h/A — 0.01. Figure [3] displays the tunneling conductance for Rd/Rb = 5 
and Eth/A = 0.01 and various h/A with (a) Z = 3 and (b) Z = 0. In this case we also find resonant ZBCP and 
splitting of the peak as in Fig. ^ The corresponding DOS of Fig. Ola) is shown in Fig. 0|for (a) h/A = and (b) 
h/A = 0.01. For h = 0, a sharp dip appears at zero energy. For finite energy the DOS is almost unity and spatially 
independent. For h/A — 0.01 a peak occurs at zero energy in the range of x near the DF/N interface. We can find 
similar structures in the corresponding conductance as shown in Fig. [31 The DOS around zero energy is strongly 
suppressed at the DF/S interface [x — L) compared to the one in Fig. [21 

Let us study the junctions with relatively large Thouless energy In this case, tunneling conductance is insensitive 
to the change of Z. In Fig. [5] we show the tunneling conductance and corresponding DOS for Z = 3, Rd/Rb = 1, 
Eth/A = 10 and various h/A. We find the broad peak of the conductance by the resonant proximity effect for 
h/A = 5 in Fig. [51(a). For h/A = 0, the DOS has a gap-like structure as shown in Fig. El(b) while for h/A = 5 it 
has a zero-energy peak as shown in Fig. [3(c). Similar plots are shown in Fig. [B] for Z = 3, Rd/Rb = 5, Eth/A = 10 
and various h/A. We find the broad ZBCP by the resonant proximity effect for h/A = 10 in Fig. [^l(a). The DOS 
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(a) h/A = 




FIG. 2: (color online) Normalized DOS for s-wave superconductors with Z = 3, Rd/Rb ~ 1 and Eth/A = 0.01. 

for h/A = has a gap- like structure as shown in Fig. [S] (b). For h/ A = 10 a zero-energy peak appears as shown in 
Fig. El (c). 

Before ending this subsection we will look at the spatial dependence of the proximity parameter, 9. FigureQdisplays 
the spatial dependence of Re# and Im# for majority spin at zero energy. We choose the same parameters as those in 
Fig0(a) and FigEl(a) for (a), (b) and (c), (d) in Fig. respectively. For the appearance of the DOS peak, large value 
of \mB is needed because the normalized DOS is given by Recos(#) = cos(Re(9)) cosh(Im(9)). When the resonant 
conditions are satisfied, Im(9 has an actually large value as shown in Fig. (b) and (d). Otherwise we can see the 
damped oscillating behavior of the proximity parameter. In contrast to Im6>, Kc9 becomes suppressed with increasing 
h/A independently of the resonant proximity effect (Fig. 0(a) and (c)). 

B. Junctions with d-wave superconductors 

In this subsection, we focus on the ci-wave junctions both for weak and strong proximity regimes. In this case, 
depending on the orientation angle a, the proximity effect is drastically changed: as a increases the proximity effect 
is suppressed^2ii&. For a = we can expect similar results to the s-wave junctions since proximity effect exists 
while the MARS is absent. On the other hand, the tunneling conductance for large a is almost independent of h/A. 
Especially, the conductance is independent of h for a/ir = 0.25 due to the complete absence of the proximity effect. 
Two different mechanisms of formation of ZBCP exist in DF/D junctions: the ZBCP caused by the resonant proximity 
effect peculiar to a ferromagnet and the ZBCP caused by the MARS located at DF/D interface. When a increases, 
MARS are formed and at the same time the proximity effect becomes weakened. Therefore the MARS provide the 
dominant contribution to the ZBCP compared to the resonant proximity effect, as will be discussed below. 
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FIG. 3: (color online) Normalized tunneling conductance for s-wave superconductors with Rd/Rb = 5 and Eth/ A = 0.01. 



First we choose the weak proximity regime where the resonant condition is h/A — 0.005. Figure [S] displays the 
tunneling conductance for Z = 3, Rd/Rb = 1 and various a with (a) Erh/A = 0.01 and h/A = 0, (b) Eth/ A = 0.01 
and h/A = 0.005, (c) E Th /A = 10 and h/A = 0, and (d) E Th /A = 10 and h/A = 5. For E Th /A = 0.01 and h = 
ZBCD appears for a/n = due to the proximity effect as in the case of the s-wave junctions while ZBCP appears for 
a/n = 0.25 due to the formation of the MARS (Fig. El (a)). For E Th /A = 0.01 and h/A = 0.005, the height of the 
ZBCP by the resonant proximity effect exceeds the one by MARS for a/ir — 0.25 (Fig. |Hl(b)). Since in the ballistic 
junctions, the ZBCP for a/ir = 0.25 is most strongly enhanced22^ii^, this ZBCP by the resonant proximity effect 
in DF is a remarkable feature. Such a feature is also expected for a larger magnitude of Eth- For Eth/ A = 10 and 
h = 0, a V-like shape of the conductance appears for a/ir — while ZBCP appears for a/ir = 0.25 (Fig. |Hl(c)). In 
this case, by choosing h/A — 5, a broad peak by the resonant proximity effect appears for a/ it = and its height 
exceeds the one for a/n = 0.25 (Fig. |U(d)). 

We also study the DOS of the DF for the same parameters as those in Fig. [S] (d) with (a) a/n = and (b) 
a/n — 0.125 in Fig. [!|] For a/n = a zero-energy peak appears as in the case of s-wave junctions. With increasing 
a the DOS around zero energy becomes suppressed due to the reduction of the proximity effect. The extreme case is 
a/n = 0.25, where the DOS is always unity since the proximity effect is completely absent. 

Next we consider the junctions in the strong proximity regime. Figure 1101 shows the tunneling conductance for 
Z = 3, R d /R b = 5 and various a with (a) E Th /A = 0.01 and h/A = 0, (b) E Th /A = 0.01 and h/A = 0.01, (c) 
E Th /A = 10 and h/A = and (d) E Th /A = 10 and h/A = 10. In this case we also find the ZBCP for a = caused 
by the resonant proximity effect. This ZBCP becomes suppressed as a increases, as shown in Figs. llOf b - ) and (d). 

The corresponding DOS of the DF for Fig. HOf d) is shown in Fig. ^2 The line shapes of the LDOS at x — are 
qualitatively similar to the tunneling conductance. The DOS at the DF/S interface (x = L) is drastically suppressed 
as compared to the one in Fig. |5J 
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FIG. 4: (color online) Normalized DOS for s-wave superconductors with Z = 3, Rd/Rb = 5 and i?Th/A = 0.01. 
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FIG. 5: (color online) Normalized tunneling conductance and corresponding DOS for s-wave superconductors with Z = 3, 
Rd/Rb = 1 and E Th /A = 10. 




FIG. 6: (color online) Normalized tunneling conductance and corresponding DOS for s-wave superconductors with Z = 3, 
Rd/Rb = 5 and E Th /A = 10. 
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FIG. 7: (color online) Spatial dependence of Re# and YmB for s-wave superconductors with Z = 3, Erh/A = 0.01. Rd/Rb = 1 
(left panels) and Rd/Rb = 5 (right panels). 
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FIG. 9: (color online) Normalized DOS for d-wave superconductors with Z — 3, Rd/Rb = 1, Eth/A. = 10 and h/A — 5. (a) 
q/7T = and (b)a/?r = 0.125. 




FIG. 10: (color online) Normalized tunneling conductance for d-wave superconductors with Z = 3 and Rd/Rb = 5. 
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FIG. 11: (color online) Normalized DOS for d-wave superconductors with Z — 3, Rd/Rb = 5, Eth/A = 10 and h/A = 10. (a) 
q/7T = and (b)a/vr = 0.125. 
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IV. CONCLUSIONS 



In the present paper, a detailed theoretical study of the tunneling conductance and the density of states in normal 
metal / diffusive ferromagnet / s- and d-wave superconductor junctions is presented. We have clarified that the 
resonant proximity effect strongly influences the tunneling conductance and the density of states. There are several 
points which have been clarified in this paper. 

1. For s-wave junctions, due to the resonant proximity effect, a sharp ZBCP appears for small Etk while a broad 
ZBCP appears for large Etk ■ We have shown that the mechanism of the ZBCP in DF/S junctions is essentially 
different from that in DN/S junctions and is due to the strong enhancement of DOS at a certain value of the exchange 
field. As a result, the magnitude of ZBCP in DF/S can exceed its normal state value in contrast to the case of DN/S 
junctions. 

2. For d- wave junctions at a — 0, similar to the s-wave case, the sharp ZBCP is formed when the resonant condition 
is satisfied. At finite misorientation angle a, the MARS contribute to the conductance when R^/Rb <C 1 and Z 3> 1. 
With the increase of a the contribution of the resonant proximity effect becomes smaller while the MARS dominate 
the conductance. As a result, for sufficiently large a ZBCP exists independently of whether the resonant condition 
is satisfied or not. In the opposite case of the weak barrier, Rd/Rb ^> 1, the contribution of MARS is negligible and 
ZBCP appears only when the resonant condition is satisfied. 

An interesting problem is a calculation of the tunneling conductance in normal metal / diffusive ferromagnet / 
p-wave superconductor junctions because interesting phenomena were predicted in diffusive normal metal / p-wave 
superconductor junctions^. We will address this problem in a separate study. 
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